Efficient implementation of high order inverse Lax-Wendroff boundary treatment for conservation laws
نویسندگان
چکیده
In [18], two of the authors developed a high order accurate numerical boundary condition procedure for hyperbolic conservation laws on a Cartesian mesh, which allows the computation using high order finite difference schemes on Cartesian meshes to solve problems in arbitrary physical domains whose boundaries do not coincide with grid lines. This procedure is based on the so-called inverse Lax-Wendroff (ILW) procedure for inflow boundary conditions and high order extrapolation for outflow boundary conditions. However, the algebra of the ILW procedure is quite heavy for two dimensional (2D) hyperbolic systems, which makes it difficult to implement the procedure for order of accuracy higher than three. In this paper, we first discuss a simplified and improved implementation for this procedure, which uses the relatively complicated ILW procedure only for the evaluation of the first order normal derivatives. Fifth order WENO type extrapolation is used for all other derivatives, regardless of the direction of the local characteristics and the smoothness of the solution. This makes the implementation of a fifth order boundary treatment practical for 2D systems with source terms. For no-penetration boundary condition of compressible inviscid flows, a further simplification is discussed, in which the evaluation of the tangential derivatives involved in the ILW procedure is avoided. We test our simplified and improved boundary treatment for Euler equations with or without source terms representing chemical reactions in detonations. The results demonstrate the designed fifth order accuracy, stability, and good performance for problems involving complicated interactions between detonation/shock waves and solid boundaries.
منابع مشابه
Stability analysis of the inverse Lax-Wendroff boundary treatment for high order upwind-biased finite difference schemes
In this paper, we consider linear stability issues for one-dimensional hyperbolic conservation laws using a class of conservative high order upwind-biased finite difference schemes, which is a prototype for the weighted essentially non-oscillatory (WENO) schemes, for initial-boundary value problems (IBVP). The inflow boundary is treated by the so-called inverse Lax-Wendroff (ILW) or simplified ...
متن کاملInverse Lax-Wendroff procedure for numerical boundary treatment of hyperbolic equations
Abstract We discuss a high order accurate numerical boundary condition for solving hyperbolic conservation laws on fixed Cartesian grids, while the physical domain can be arbitrarily shaped and moving. Compared with body-fitted meshes, the biggest advantage of Cartesian grids is that the grid generation is trivial. The challenge is however that the physical boundary does not usually coincide wi...
متن کاملInverse Lax-Wendroff procedure for numerical boundary conditions of conservation laws
We develop a high order finite difference numerical boundary condition for solving hyperbolic conservation laws on a Cartesian mesh. The challenge results from the wide stencil of the interior high order scheme and the fact that the boundary intersects the grids in an arbitrary fashion. Our method is based on an inverse Lax-Wendroff procedure for the inflow boundary conditions. We repeatedly us...
متن کاملAn inverse Lax-Wendroff method for boundary conditions applied to Boltzmann type models
In this paper we present a new algorithm based on a Cartesian mesh for the numerical approximation of kinetic models on complex geometry boundary. Due to the high dimensional property, numerical algorithms based on unstructured meshes for a complex geometry are not appropriate. Here we propose to develop an inverse Lax-Wendroff procedure, which was recently introduced for conservation laws [21]...
متن کاملA high order moving boundary treatment for compressible inviscid flows
We develop a high order numerical boundary condition for compressible inviscid flows involving complex moving geometries. It is based on finite difference methods on fixed Cartesian meshes which pose a challenge that the moving boundaries intersect the grid lines in an arbitrary fashion. Our method is an extension of the so-called inverse Lax-Wendroff procedure proposed in [16] for conservation...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- J. Comput. Physics
دوره 231 شماره
صفحات -
تاریخ انتشار 2012